C::: ■ FIRST RUN II MEASUREMENT OF THE W BOSON MASS WITH CDF 
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, The CDF collaboration has analyzed ~200 pb~^ of Tevatron Run 11 data taken between 

\^ • February 2002 and September 2003 to measure the W boson mass. With a sample of 63964 

' W ^ eu decays and 51128 W —* ^lv decays, we measure Mw = 80413±34(stat)±34(syst) 

, MeV/c^. The total measurement uncertainty of 48 MeV/c^ makes this result the most precise 

' single measurement of the W boson mass to date. 




1 Introduction 

The W boson mass is an important Standard Model (SM) parameter. It receives self-energy 
corrections due to vacuum fluctuations involving virtual particles. Thus the W boson mass 
probes the particle spectrum in nature, including particles that have yet to be observed directly. 
The hypothetical particle of most immediate interest is the Higgs boson. The interaction of the 
other SM particles, in particular the W and Z gauge bosons, with the Higgs field is thought to 
impart mass to the SM particles. Thus the Higgs boson plays a critical role in the SM and it 
is very interesting to obtain a Higgs mass prediction. The W boson mass can be calculated at 
tree level using the three precise measurements of the Z boson mass, the Fermi coupling Gp 
and the electromagnetic coupling Oem- In order to extract information on new particles, we 
need to account for the radiative corrections to My/ due to the dominant top-bottom quark 
loop diagrams. For fixed values of other inputs, the current uncertainty on the top quark mass 
measurement 170.9ibl.8 GeV/c'^ — corresponds to an uncertainty in its W boson mass correction 
of 11 MeV/c^. Measurements of the W boson mass from Run I of the Tevatron and LEP with 
uncertainties of 59 MeV/c^^^and 33 MeV/c^^ respectively, yield a world average of 80392zb29 
MeV/c^^. It is clearly profitable to reduce the W boson mass uncertainty further as a means 
of constraining the Higgs boson mass. 



2 Measurement Strategy 



At the Tevatron, W bosons are mainly produced by valance quark-antiquark annihilation, with 
initial state gluon radiation (ISR) generating a typical transverse boost of 0(10 GeV). The 
transverse momentum (pip) distribution of the decay lepton has a characteristic Jacobian edge 
whose location, while sensitive to the W boson mass, is smeared by the transverse boost of the 
W boson. The neutrino transverse momentum (pj,) can be inferred by imposing px balance 

in the event. The transverse mass, defined as rriT = \j2pipPj,{l — cos[(j)^ — (f)^]), includes both 
measurable quantities in the W decay and provides the most precise quantity to measure the W 
boson mass. We use the mr, Pt and pj^ distributions to extract Mw. These distributions do not 
lend themselves to analytic parameterizations, which leads us to use a Monte Carlo simulation 
to predict their shape as a function of M\y. These lineshape predictions depend on a number of 
physical and detector effects, which we constrain from control samples or calculation. Important 
physical effects include internal QED radiation, the intrinsic W boson transverse momentum, 
and the proton parton distribution functions. Detector effects include external bremsstrahlung 
and ionization energy loss in the detector material, tracker momentum scale, calorimeter energy 
scale, resolution of the tracker and calorimeter, and the detector fiducial acceptance. In order 
to model and study these effects at the level of a part in 10^, we have developed a sophisticated, 
first-principles fast Monte Carlo simulation of the CDF detector. The W boson mass is extracted 
by performing a binned maximum-likelihood fit to the rriT, Pt and p^ distributions. We generate 
800 templates as a function of Mw between 80 GeV/c^ and 81 GeV/c^. 

3 Momentum and Energy Scale Calibration 

The key aspect of the measurement is the calibration of the lepton momentum. The charged 
lepton momentum is measured in a cylindrical drift chamber called the Central Outer Tracker 
(COT). The electron energy is measured using the central electromagnetic (EM) calorimeter and 
its angle measurement is provided by the COT trajectory. The COT track measurement sets the 
momentum scale for this analysis. The internal alignment of the COT is performed using high- 
Pt cosmic rays that traverse diametrically the entire drift chamber. The momentum scale is set 
by measuring the J/^ and T(15') masses using the dimuon mass peaks. The J/^ sample spans 
a range of muon pj- (2-10 GeV/c), which allows us to tune our ionization energy loss model such 
that the measured mass is independent of muon pT- We obtain consistent calibrations from the 
J/^, T(IS') mass fits shown in Fig. [1] (left). The momentum scale extracted from the Z fi^ 




Figure 1: Left: Momentum scale summary: Ap/p vs l/pr for J/^', T(IS') and Z boson dimuon samples. The 
dotted line represents the independent uncertainty between J/^' and T(IS'). Error bars are statistical only. 
Right: Energy scale calibration using E/p distribution from W eu events. 



mass fit, shown in the same figure, is also consistent, albeit with a larger, statistics-dominated 
uncertainty. The tracker resolution is tuned on the observed width of the T(15) and Z boson 
mass peaks. 

Given the tracker momentum calibration, we fit the peak of the Ejp distribution of the signal 
electrons in the W ^ eu sample (Fig. [1] right) in order to calibrate the energy measurement 
of the electromagnetic (EM) calorimeter. The energy scale is adjusted such that the fit to 
the peak returns unity. The model for radiative energy loss is constrained, by comparing the 
number of events in the radiative tail of the E/p distribution. The calorimeter energy calibration 
is performed in bins of electron Et to constrain the calorimeter non-linearity. The calibration 
yields a Z ^ ee mass measurement of Mz = 91190ib67stat MeV/c^, in very good agreement with 
the world average (91187.6 ±2.1 MeV/c^E!); we obtain the most precise calorimeter calibration 
by combining the results from the E/p method and the Z ^ ee mass measurement. The EM 
calorimeter resolution model is tuned on the widths of the E/p peak and the Z — > ee mass peak, 
separately for non-radiative and radiative electrons. 

4 Hadronic Recoil Calibration 

All particles recoiling against the W or Z boson are collectively referred to as the recoil. The 
recoil is computed as the vector sum of transverse energy over all calorimeter towers, where the 
towers associated with the leptons are explicitly removed from the calculation. The response of 
the calorimeter to the recoil is described by a response function which scales the true recoil mag- 
nitude to simulate the measured magnitude. The hadronic resolution receives contributions from 
ISR jets and the underlying event. The latter is independent of the boson transverse momentum 
and modeled using minimum bias data. The recoil response and resolution parameterizations 
are tuned on the mean and rms of the pT^-imbalance in Z ^ II events as a function of boson 
Pt- We define the r/ axis to be the geometric bisector of the two leptons and the ^ axis to be 
perpendicular to rj. We project the vector pj^-balance onto the r] and ^ axes and compare the 
data distribution to the simulation. Fig. [2] (left) shows the mean of the pT-balancing in Z — > ee 
events as a function of Z boson pT- The quantity p^ is computed from the EM clusters of the 
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Figure 2: Left; Mean of the pr-balancing as a function of pt in Z ee events. Right: Projection of recoil along 

the lepton direction in W ^ j-if events. 

two Z boson decay electrons and is computed from the recoil vector in the calorimeter. We 
cross-check the recoil model using W and Z boson data which show good agreement and validate 
the model. A very sensitive quantity to cross-check the recoil model in W boson events is the 
projection of the recoil along the lepton direction (ti||). Fig. [5] (right) shows U|| for W ^ pv 
events. We find good agreement between the data and the simulation. 



5 Event Generation 



We generate W and Z events with resbosSI, which captures the QCD physics and models the W 
Pt spectrum, resbos calculates the quintuple differential production cross section dqTd!^Ct ^ 
where Q, y and qt are the boson invariant mass, rapidity, and transverse momentum respectively, 
and is the solid angle element in the decay lepton direction. The resbos parametrization of 
the non-pertubative form factor is tuned on the dilepton pT distribution in the Z boson sample. 
Photons radiated off the final-state leptons (FSR) are generated according to wgradEI. The 
FSR photon energies are increased by 10% (with an absolute uncertainty of 5%) to account 
for additional energy loss due to two-photon radiation! WGRAD is also used to estimate the 
uncertainty due to QED radiation from initial state (ISR) and interference between ISR and 
FSR. We use the CTEQ6M^ set of parton distribution functions at NLO and apply their 
uncertainties to evaluate the systematic uncertainty on the W boson mass. 



6 Backgrounds 

Backgrounds passing the event selection have different kinematic distributions from the W signal 
and are included in the template fit according to their normalizations. Backgrounds arise in the 
W boson samples from misidentified jets containing high-p^^ tracks and EM clusters, Z ^ II 
where one of the leptons is not reconstructed and mimics a neutrino, W — > rz^, kaon and pion 
decays in flight (DIF), and cosmic ray muons. The latter two are backgrounds in the muon 
channel only. Jet, DIF, and cosmic ray backgrounds are estimated from the data to be together 
less than 0.5%. The W ^ tv background is 0.9% for both channels, and the Z ^ II \s 6.6% 
(0.24%) in the muon (electron) channel, as estimated from Monte Carlo samples generated with 
PYTHIA and a detailed GEANT-based detector simulation. 



7 Results and Conclusions 

The fits to the three kinematic distributions my, Pt and p^ in the electron and muon channels 
give the W boson mass results shown in Table [TJ 

Table 1: Fit results from the distributions used to extract Mw with uncertainties. 



Distribution 


W boson mass (MeV/c^) 


xVdof 




80493±48,tat±39,y,t 


86/48 


Me) 


80451±58,tat±45,j,,t 


63/62 


pUe) 


80m±57 stat±54,yst 


63/62 


mTilJ-,1^) 


80M9±5Astat±27syst 


59/48 




80321±66,tat±40,y,t 


72/62 


pW 


80396±66,tat±46,y,j 


44/62 



The transverse mass fit for the W ^ ev channel is shown in Fig. [3] (left) and for the 
W ^ channel in Fig. [3] (right). The uncertainties for the rriT fits in both channels are 
summarized in Table [2j We combine the six W boson mass fits including all correlations to 
obtain Mvi/=80413±34(stat)±34(syst) MeV/c^. With a total uncertainty of 48 MeV/c^, this 
measurement is the most precise single measurement to date. Inclusion of this result increases 
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Figure 3: Left: Transverse mass fit in the electron decay channel. Right: Transverse mass fit in the muon decay 
channel. The histograms at the bottom of the distributions represent the backgrounds. 



the world average W boson mass to M(4'=80398±25 MeV/c^^ reducing its uncertainty by 15%. 
The updated world average impacts the global precision electroweak fits, reducing the preferred 
Higgs boson mass fit by 6 GeV/c^ to M//=76]'l24 GeV/c^. The resulting 95% CL upper limit on 
the Higgs mass is 144 GeV/c2 (182 GeV/c^) with the LEP II direct limit includedjexcluded)'^ 
The direction of this change has interesting theoretical implications: as FigUpl shows, the 
M]Y vs Mtop ellipse moves a little deeper into the light-Higgs region excluded by LEP II and into 
the region favored by the minimal supersymmetry model (MSSM). While this is a one-sigma 
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Figure 4: Constraint on the Higgs boson mass from direct Mw and Mtop measurements at 68%, 95% and 99.7% 

CL, along with SM and MSSM calculations. 

effect, it arouses further interest in higher precision measurements of Mw (and Mtop)- Most 
of the systematic uncertainties in this measurement (Table [2]) are limited by the statistics of 
the calibration samples used. Further improvements in the detector model and the production 
and decay model (e.g. QED radiative corrections) are likely to shrink the other systematic 



Table 2: Systematic and total uncertainties for the transverse mass fits. The last column shows the correlated 
uncertainties, the last row is the quadrature sum of statistical and systematic uncertainty. 



Systematic (MeV/c^) 


W ^eu 


W ^ fxu 


Common 




30 


17 


17 


Lepton Energy Resolution 


9 


3 





Recoil Energy Scale 


9 


9 


9 


Recoil Energy Resolution 


7 


7 


7 


Selection Bias 


3 


1 





Lepton Removal 


8 


5 


5 


R a rl<" p" r o 1 Tn H s 


8 


9 





Pt{W) Model 


3 


3 


3 


Parton Distributions 


11 


11 


11 


QED radiation 


11 


12 


11 


Total Systematics 


39 


27 


26 


Total Uncertainty 


62 


60 


26 



uncertainties as well. CDF has now accumulated an integrated luminosity of about 2 fb^ and 
we look forward to a boson mass measurement with precision better than the current world 
average of 25 MeV/c^, with the dataset already in hand. 
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